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In simulations of a lOPW laser striking a solid we demonstrate the possibility of producing a 
pure electron-positron plasma by the same processes as those thought to operate in high-energy 
astrophysical environments. A maximum positron density of lO^^'m"'^ can be achieved, seven orders 
of magnitude greater than achieved in previous experiments. Additionally, 35% of the laser energy 
is converted to a burst of gamma-rays of intensity lO^^Wcm"^, potentially the most intense gamma- 
ray source available in the laboratory. This absorption results in a strong feedback between both 
pair and 7-ray production and classical plasma physics in the new 'QED-plasma' regime. 



Electron-positron (e^e"*") plasmas are a prominent fea- 
ture of the winds from pulsars and black holes [H 
They result from the presence of electromagnetic fields 
strong enough to cause non-linear quantum electrody- 
namics (QED) reactions [3] in these environments lead- 
ing to a cascade of e~e+ pair production [1]. These fields 
can be much lower than the Schwinger field for vacuum 
breakdown [S] if they interact with highly relativistic elec- 
trons (7 >> 1) |3]. Non-linear QED has been probed 
experimentally with lasers in two complementary ways: 
(1) with a particle accelerator accelerating electrons to 
the necessary 7 and a laser supplying the fields [5HH]; 
or (2) with a laser accelerating the electrons and gold- 
nuclei supplying the fields [QUIT] . An alternative configu- 
ration, using next-generation high-intensity lasers to pro- 
vide both the acceleration and the fields [T^ , has the po- 
tential to generate dense e~e+ plasmas. Analytical calcu- 
lations and simulations exploring this configuration have 
shown that an overdense e~e^ plasma can be generated 
from a single electron by counter-propagating lOOPW 
lasers p!^15) . Here we will show that such a plasma 
can be generated with an order of magnitude less laser 
power by firing the laser at a solid target, putting such 
experiments in reach of next-generation lOPW lasers |16j . 

The dominant non-linear QED effects in lOPW laser- 
plasma interactions are: synchrotron gamma-ray pho- 
ton (7^) emission from electrons in the laser's electro- 
magnetic fields; and pair-production by the multipho- 
ton Breit- Wheeler process, + ^^7; — e~ -I- e+, where 
7i is a laser photon [31 [T71 [TS]. Each reaction is a 
strongly multiphoton process, the former process being 
non-linear Compton scattering, e~ -f m^yi e.^ + 
[TOl [5D], in the limit m — >■ 00. Therefore, these reac- 
tions only become important at the ultra-high intensities 
reached in lOPW laser-plasma interactions. The impor- 
tance of synchrotron emission is determined by the pa- 
rameter T]. This depends on the ratio of the electric and 
magnetic fields in the plasma to the Schwinger field [5] 
{Es = 1.3 X lO^^Vm"^). For ultra-relativistic particles 



ri = (7/£:s)|E_L + /? X cB| HZl [H]. 7 is the Lorentz 
factor of the emitting electron or positron, (3 is the corre- 
sponding velocity normalised to c and E_l is the electric 
field perpendicular to its motion. As rj approaches unity 
each emitted photon takes a large fraction (« 0.44) of 
the emitting electron's energy, and the mean free path of 
these photons to pair production is of the order of the 
laser wavelength so that many pairs are produced |12j . 
For a lOPW laser operating at an intensity of lO^'^Wcm"^ 
|E| w lO^^Vm"^. On interacting with a plasma such a 
laser accelerates electrons to a 7 of the order of several 
hundreds and so 77 approaches one. However, the geome- 
try of the interaction is crucial; for a single intense laser 
beam striking a single electron, the electron is rapidly 
accelerated to w c in the direction of propagation of the 
laser pulse. In this case E_l is almost exactly cancelled by 
V X B, 77 is reduced and pair production is dramatically 
curtailed. By contrast, in an overdense plasma the wave 
becomes evanescent and the terms do not cancel. There- 
fore, laser-solid interactions offer an attractive route to 
generating electron-positron plasmas. 

In this Letter we will present the first simulations 
of lOPW laser-solid interactions to include the relevant 
QED processes. We show that such interactions are the 
most effective way to produce e~e"'" pairs with next- 
generation lasers and that the laser is absorbed into an 
ultra-intense burst of 7-rays with high efficiency (35%). 
In order to understand these interactions it is crucial to 
resolve the complex feedback between QED and collec- 
tive plasma physics effects. Therefore, in contrast to the 
schemes described in (1) & (2) above, we describe a new 
regime in which QED processes and plasma physics are 
inseparable, which we term a 'QED-plasma'. 

In order to simulate QED-plasmas we have included 
synchrotron emission of high energy 7-ray photons and 
Breit- Wheeler pair production in the particle-in-cell 
(PIC) code EPOCH |21]- As 77 approaches unity the 
high energy of the emitted photons means that radia- 
tion must be considered discontinuously. The electrons 
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Figure 1: Pair-production by a laser of intensity 4 x 
lO^'^Wcm"'^ striking an aluminium target (snapshots at the 
end of the 30fs laser pulse). The laser (red contours) bores 
a hole into the solid target (blue). Gamma-rays (blue) and 
positrons (red dots) are generated in this interaction (inset - 
on the same scale). 
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Figure 2: Electron & positron energy spectra in the inter- 
action of a laser of intensity 4 x lO^'^Wcm"^ with solid alu- 
minium (spatially and temporally integrated). As usual the 
electron spectrum has a pronounced tail of 'fast electrons'. 
Corresponding gamma-ray spectrum (inset). 



and positrons obey the Lorentz force equation, following 
the classical worldlines as computed by the PIC code, 
until a discrete photon is emitted [52]. The recoil in 
such an event provides a discontinuous radiation reac- 
tion force [5D]. As discussed below, the discontinuous 
radiation model consists of random sampling of the syn- 
chrotron spectrum and so tends to the continuous-loss 
model [TT l |23H26| as hojh « 7meC^ (/iw/i is the energy 
of the emitted photon), i.e. as the sampling frequency 
oo. It has recently been shown that in lOPW laser- 
plasma interactions the discontinuous model yields an 
order of magnitude more e~e'^ pairs [27j . This is due 
to some electrons reaching higher energies and emitting 
a higher energy photon than the same electron experi- 
encing a continuous radiation drag force, the so-called 
'straggling' effect [22j. 

The QED processes are simulated using a Monte-Carlo 
algorithm [27] . The time at which emission events oc- 
cur is computed as follows. Each particle is assigned 
an optical depth at which it emits (r) according to 
P = 1 — e""^, where P G [0, 1] is chosen at random to 
capture the quantum fluctuations in the emission pro- 
cesses and so the straggling. The rates of photon and pair 
production, dr^/dt = {V^afcrj)/ (Act) Iq^"^ dxF{T], x)/x 
and dT±/dt = {2Trafc/ Xc)imeC^ /hujh)xT±{x)^ a-i'e then 
solved until these optical depths are reached, at which 
point the emission event occurs [27j . Here, a/ is the 
fine structure constant and Ac is the Compton wave- 
length; X = {hujii/2mec'^)\E,± + k X cB|, where _L signi- 
fies the field component perpendicular to the unit vec- 
tor in the photon's direction of motion k. Photons are 
generated with a random energy weighted by the syn- 
chrotron function F{rj,x) including Klein-Nishina correc- 
tions (18j . X controls pair production via the function 



r±(x) « Q.lQKli^{2/'ix)/x- The generated pairs are 
treated on an equivalent footing to the original electrons 
in the PIC code, the photons are treated as massless, 
chargeless macroparticles which propagate ballistically. 
The pairs are included when the PIC code calculates the 
charge and current densities on the computational grid, 
so contribute to the electromagnetic fields that are used 
to calculate the QED rates at the next timestep; ensuring 
a self-consistent simulation. 

We have performed two-dimensional EPOCH simula- 
tions of a lOPW laser striking an aluminium foil includ- 
ing the QED processes. The results are shown in Fig. 
[T] The aluminium target is 1/im thick, has a density of 
2700kgm^'^ and is assumed to be fully ionised. The target 
is represented by 1000 pseudo-electrons and 32 pseudo- 
ions per cell, with a spatial resolution of lOnm. The laser 
has wavelength A; = l//m and is linearly p-polarised. The 
pulse has an energy of 377J and a duration of 30fs, with a 
square temporal profile. It is focused to a spot of radius 
1/zm with intensity / = 4x lO^'^Wcm"^. For this laser in- 
tensity the electron density of fully ionised aluminium rig 
is higher than the relativistically-corrected critical den- 
sity Tic = ^me^Qi^f / (w/ = 27rc/A/) and the plasma is 
overdense. Therefore the laser beam is refiected and the 
light-pressure of the beam bores a hole into the target as 
shown in Fig. [T] Also shown is prolific gamma-ray and 
positron production at the hole-boring front, where the 
laser is reflected. The total number of pairs produced is 
N± = 8 X 10^ (each red dot is a macroparticle represent- 
ing 2 X 10^ positrons). 

Pairs are overall electrically neutral and so readily es- 
cape the target. Thin sheets of pure electron-positron 
plasma form behind the target with a positron number 
density of lO^^m"^. An e~e^ plasma is also trapped in- 
side the hole-boring cavity with density lO^^m""^ over one 
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Figure 3: Gamma-ray intensity averaged over one laser pe- 
riod for a laser intensity of 8 x 10^"^Wcm~^ 25fs after the end 
of the incident laser pulse (when all the photons leave the 
target). 



cubic micron, forming a self-contained 'micro-laboratory' 
potentially useful for the study of such a plasma. For the 
l^m thick target the laser just breaks through the target 
at the end of the BOfs laser pulse, releasing the trapped 
pairs for probing. When the laser breaks through the 
situation reverts to that of a single electron in a sin- 
gle beam, pair production ceases and further laser en- 
ergy is wasted. The positron density is seven orders 
of magnitude higher than produced by the gold-target 
scheme described above and is high enough that collec- 
tive effects could be studied with a CO2 laser. Fig. [2] 
shows that the average positron energy of 250MeV is 
much higher than the energy of photons from which they 
originate. This suggests that the positrons are acceler- 
ated to high energy by the laser. In this case we ex- 
pect the average Lorentz factor of the positrons to be 
(7) = a^°'-K$ = 2a'°' = 2eEffymeCLUi saSOOMeV, which 
is consistent with the simulations. Ejfg is the value of 
the electric field inside the solid (as discussed below). $ 
is the sheath potential generated by fast electrons as they 
leave the target. The sheath field acts to confine the fast 
electrons (the majority species compared to positrons) 
inside the target and so accelerates positrons [H], do- 
ing work approximately equal to the fast electron energy 
P5]). In practise lasers can have a long timescale pre- 
pulse of lower intensity than the main pulse. Such a pre- 
pulse may heat the target and cause it to expand prior 
to the arrival of the main pulse, generating a pre-plasma. 
Additional simulations, similar to that discussed above, 
show that a small (e-folding distance = l/xm) pre-plasma 
does not dramatically reduce the number of gamma-ray 
photons and pairs generated. In fact, for the parameters 
explored here a pre-plasma actually enhances gamma- 
ray production by 10%. A full exploration of this en- 
hancement over all parameter space, as well as the role 
of pre-plasma in gamma-ray and pair production for a 
laser pulse at oblique incidence is beyond the scope of 
this paper. 



Synchrotron gamma-ray photons are generated prolif- 
ically in the laser-solid interaction [29]. At a laser inten- 
sity of 8 X lO^'^Wcm"^, a burst of gamma-rays of average 
intensity 8 x lO^'^Wcm"^ is produced at the rear of a 1/im 
thick Al target. This is shown in Fig. [sj 10^^ gamma-ray 
photons with an average energy of 16McV are produced. 
The conversion fraction of laser to gamma-ray energy 
is 0.35. Note that synchrotron emission dominates over 
the more usual bremsstrahlung emission (not included in 
the simulation) during the laser pulse. The synchrotron 
emission occurs during the pulse duration, whereas the 
bremsstrahlung cooling time (several picoseconds) is sub- 
stantially longer. In the simulation the synchrotron 
gamma-ray emission is contained within a cone half-angle 
of (f>sim — 80°. which is consistent with the relativisti- 
cally boosted angle (j)boost = cos~^ {vhb/c). Here vhb 
is the speed of the hole-boring front (discussed further 
below). The substantial energy loss to gamma-ray emis- 
sion profoundly alters the energy budget of the laser-solid 
interaction and so the plasma physics processes. The av- 
erage electron energy is reduced from 41MeV to 21MeV 
|30) . The average ion energy is not strongly modified but 
the spectrum is substantially modified by the synchrotron 
emission; altering it from a single peak at 3.5GeV (the 
hole boring model, described below, predicts 3GeV) to 
two peaks at 1.5GeV & 4.5GeV. The reduced peak is due 
to the reduction in the reflected laser intensity, reducing 
the strength of the pistoning of the ion surface from 21 /c 
for perfect reflection to I/c for perfect absorption. The 
enhanced energy peak occurs as the laser breaks through, 
as show by Tamburini et al & Chen et al [301 131] • 

We have shown that gamma-ray emission alters the 
energy budget of laser-solid interactions and so the clas- 
sical plasma physics. Conversely, the rates of the QED 
reactions are strongly modified by the plasma physics 
processes, closing the feedback loop which is the defin- 
ing feature of QED-plasmas. The modification of the 
QED rates can be estimated by employing the analytical 
model of Bell & Kirk [T^. Here the controlling param- 
eter 77 is expressed implicitly in terms of 124, the laser 
intensity in units of lO^^Wcm"^, and A^m, the laser 
wavelength in microns as I24 — 2.75?]'' -I- 0.28?7/Aprn- In 
laser-solid interactions three plasma effects reduce i] and 
consequently the QED rates. (1) Relativistic hole-boring 
[32]: when / < 8 x lO^^Wcm"^ the laser reflects from the 
overdense solid's hole-boring surface, which is moving at 
relativistic speed vhb- Where vhb/c = •\/S/(l -I- -s/S) 
and S = 1 1 pt? is the dimensionless pistoning parame- 
ter. The energy of the accelerated ions (of mass m^) is 
2SmiC^/(l -I- 2-\/S)- In the rest-frame of the hole bor- 
ing surface, the intensity and wavelength are modified 
by the relativistic Doppler effect to I24,hb & A^m,j?B 
(2) The skin-effect: the maximum value of the electric 
field in the evanescent wave inside the solid is reduced 
to i;""' = 2(nc/nejjB)i/2£:™^^. is the peak laser 

electric field and rieHB is the electron number density in 
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the hole-boring frame. This can be included by mod- 
ifying I24.HB to lit^HB = hi.HBnc/rieHB- Maximum 
I24HB achieved by reducing the target density such 
that it is just above the relativistic critical density at 
the incident laser intensity. (3) Self-induced transparency 
[33] : for / > 8 X lO^^Wcm^^ the solid target begins to 
become transparent, the situation approaches that of a 
laser interacting with a single electron and the rate of 
pair production is strongly reduced. 

The equation for rj including plasma physics effects 
«b) is IltHB = 2.75(r;-^)4 + 0.28^jf's/\^mMB. We 
can solve this numerically to obtain ?7|f^ in a given 
laser-solid interaction and so estimate the number of 
gamma-ray photons produced and their energy. To do 
this we use the rate equation for dr^/dt given above 
with 77 ^ r]jf^s and F{ri,x) ~^ fmono{4:X/^V^ = 2/) = 
{8TT/9V3)S{y — 0.29). The latter corresponds to assum- 
ing that the emitted photons are monochromatic with 
energy {hvh) = 0.44r/|^°jj(7)mec2 [ig. The number of 
gamma-ray photons produced per electron per laser pe- 
riod is then given by A^^ = 6.42a/ (7) [T2]. (7) is the aver- 
age Lorentz factor of the electrons and can be estimated 
by (7) ~ a'^°^- For a laser of intensity 8 x lO^^Wcm"^ 
focused onto a solid aluminium target jylfg « 0.4 and so 
« 4x 10^'^ and huh ~ 25MeV. These are in reasonable 
agreement with the simulation results presented above. 

An alternative configuration for pair production, re- 
cently investigated by Nerush et al [TS] & Elkina et al 
|34j . is the interaction of counter-propagating lasers in 
an underdense gas. In this case plasma physics effects 
do not reduce the QED rates, but the plasma density is 
much lower. In order to compare these configurations 
we performed one-dimensional EPOCH simulations of: 

(1) a laser of intensity / striking a solid, semi- infinite 
(to avoid complicating break-through effects), Al target; 

(2) counter-propagating lasers of intensity 1/2 in an un- 
derdense hydrogen gas-jet. For / < 8 x lO^^Wcm^^ 
more pairs are produced by the solid target configura- 
tion (above this intensity the aluminium target becomes 
transparent). The 10^ times denser plasma outweighs 
the 10^ — 10''' times reduced rate of reaction for the solid. 
This rate reduction can be estimated analytically. The 
number of pairs produced per electron per laser period 
is N-y{l — er^'^'^) [12], where (r) is the photon optical 
depth for absorption over a distance A;. Here (t) — 
12.8/24e-4/3(x> and (x) = {{hwh) /2me^){EI^\,l E^). For 
a / = 8 X lO^^Wcm^^ laser-aluminium interaction the re- 
duction in rj leads to a reduction in N± by a factor of 10^, 
in good agreement with the simulations. For intensities 
of the order of 10^'*Wcm~^ (expected to be reached by 
lOOPW class lasers) the gas-jet configuration produces 
more pairs. In contrast to the solid a large fraction of 
the pairs generated go on to produce additional pairs, the 
reaction runs away and a cascade of antimatter produc- 
tion ensues. This is in good agreement with the results 



of Nerush et al [T5] . 

In conclusion, we have shown that lOPW laser-solid 
interactions will generate dense electron-positron plas- 
mas and ultra-intense bursts of gamma-rays, relevant to 
the laboratory study of pair production in high-energy 
astrophysical environments. In contrast to the other 
laser-based positron production schemes mentioned, we 
have shown that for lOPW laser-solid interactions there 
is a strong feedback between QED processes and plasma 
physics, leading to the new regime of QED-plasma 
physics. An understanding of future experiments in 
this regime will be impossible without a self-consistent 
model including the interplay between QED and classi- 
cal plasma physics as discussed here. 
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